Slow dynamics and aging of a confined granular flow 
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We present experimental results on slow flow properties of a granular assembly confined in a 
vertical column and driven upwards at a constant velocity V. For monodisperse assemblies this 
study evidences at low velocities (f <V<fOO fim/s) a stiffening behaviour i.e. the stress necessary 
to obtain a steady sate velocity increases roughly logarithmically with velocity. On the other hand, 
at very low driving velocity (V < 1 fim/s), we evidence a discontinuous and hysteretic transition to a 
stick-slip regime characterized by a strong divergence of the maximal blockage force when the velocity 
goes to zero. We show that all this phenomenology is strongly influenced by surrounding humidity. 
We also present a tentative to establish a link between the granular rheology and the solid friction 
forces between the wall and the grains. We base our discussions on a simple theoretical model and 
independent grain/wall tribology measurements. We also use finite elements numerical simulations 
to confront experimental results to isotropic elasticity. A second system made of polydisperse 
assemblies of glass beads is investigated. We emphasize the onset of a new dynamical behavior, i.e. 
the large distribution of blockage forces evidenced in the stick-slip regime. 

PACS numbers: 45.70.-n,46.55.+d,81.05.-Rm 
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I. INTRODUCTION 



Granular flows are presently at the focus of many at- 
tentions Q] . The classical model for granular media sta- 
bility was proposed at the end of the XVIIIth century by 
Ch. A. de Coulomb who revealed a strong analogy be- 
tween the failure properties of a granular assembly and 
the phenomenology of solid on solid friction. Modern de- 
velopments have elaborated sophisticated empirical ap- 
proaches around this fundamental idea 0] but so far there 
is no deep physical understanding nor rigorous derivation 
describing the passage from the granular level descrip- 
tion to a set of evolution equations involving macroscopic 
quantities such as stress, strain or packing fraction. Fur- 
thermore, there is an additional difficulty to apprehend 
complex behaviors such as aging under stress or to ac- 
count for the influence of external parameters such as 
surrounding humidity which effects are often noticed in 
practice. Note that similar questions are still under ac- 
tive consideration in the field of tribology 0, 0, H , but 
in the case of granular assemblies, a supplementary dif- 
ficulty lies in the fragile character of granular structures 
which can be easily modified under the action of external 
constraints. 

Dynamical behavior of slowly driven granular materi- 
als was investigated by different groups both in compres- 
sion and/or in shearing experiments 0, ■ Here we 
present an experimental situation of the same type, but 
in a quite different geometry. We investigate the rheology 
of a granular assembly confined in a cylindrical column 
and pushed vertically from the bottom. The resistance to 
vertical motion as well as the blocking/unblocking tran- 
sitions reveals a phenomenology possibly shared by many 
confined granular assemblies. Note that this column con- 



figuration may help to understand several practical sit- 
uations like pipe flows ^(j> compaction under stress or 
dense granular paste extrusion. A previous investigation 
of the same display was made in 2D an d also in 3D 
as a preliminary report, and more recently see refs 
|l3l fl4| . Those contributions have shown a rich phe- 
nomenology partly sorted by the solid friction properties 
of the grains and the boundaries. The importance of 
surrounding humidity was also evidenced |l2^ as it would 
strongly influence the rheology of the granular column. 
In this report, we push further the investigation as we 
change the column material, the beads characteristics 
and size dispersity, under various humidity conditions. 
We also propose to compare our data with the outcome 
of a simple numerical model of isotropic elasticity. 

In the classical situation of newtonian fluid pushed by 
a piston (typically in a syringe), one would obtain a re- 
lation between pressure and flow rate, which is charac- 
terized by a fluid constitutive parameter: the viscosity. 
For an isotropic elastic medium the resistance to push- 
ing would depend on the material Poisson's ratio and wall 
frictional properties 0, 0] (we will detail this question 
further in the text). For a granular material, the situa- 
tion is a priori more complex since the piling structure 
can be modified so as to adapt to the external constraints. 
Friction at the walls may also create internal granular re- 
circulation flows Therefore it is an important but 
difficult matter to clarify the rheology of this system by 
sorting the respective influence of grain-boundary fric- 
tion versus bulk structural changes. Along this line we 
started with a simple situation of weakly frictional steel 
beads with frictional boundaries. A second situation is 
studied where the granular material is made of rugous 
polydisperse glass beads. 
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FIG. 1: Sketch of the experimental display. Inset: sketch of 
the slider. 



II. DISPLAY 

Most experiments are performed on dry, non cohesive 
and monodisperse steel beads of diameter d — 1.58 mm 
piled into a vertical duralumin or brass cylinder of diame- 
ter D = 36 mm. We also perform a series of experiments 
on a miscellany of glass beads of three diameters (1.5 
mm, 2 mm and 3 mm) with equal volume of each kind, 
in a PMMA cylinder. The column is closed at the bot- 
tom by a movable brass piston avoiding contact with the 
column (diameter mismatch is 0.5 mm). A force probe 
of stiffness k = 40000 N.m -1 is located under the piston 
and is pushed at a constant driving velocity V (between 
5 nm.s -1 and 100 /im.s -1 ) via a stepping motor (see Fig. 

The resistance force F encountered by the piston is 
measured as a function of time. We also monitor the rel- 
ative humidity (\) and the surrounding temperature. We 
work in the range 35% < x < 75%, as well as in dry air 
(x < 3%) and humid air (x = 90%). Actually, except for 
the dry and humid situation, we do not regulate this last 
parameter (x) but we record its values close to the ex- 
perimental set-up. We obtain dry air by having a weak 
air flux (250 ml.min -1 ) flow into a cylinder filled with 
silica. We obtain humid air by making a weak air flux 
(250 ml.min -1 ) bubble through water. Then, for the du- 
ration of an experiment, weak fluxes of humid or dry air 
are set to flow through the column from top to bottom, 
maintaining a constant and homogeneous level of humid- 
ity. Note that a stationary relative humidity is reached in 
a few minutes. We start our experiments when this sta- 
tionary level is attained. We verified that the weak air 
flux does not perturb the system: we notice no evolution 
of the force when the flux is stopped as long as relative 
humidity remains unchanged. Temperature is kept at 
(20 ± 1) ° C. We actually find no correlation between the 
force fluctuations and the temperature variations in this 
range. 

Note that the brass and duralumin cylinders have very 
different surface properties (Fig. |2J). The duralumin 
cylinder is rough cast, and its roughness is 400 nm. The 
brass cylinder was machine-turned, and has mean rough- 
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FIG. 2: Topographies of internal surfaces of duralumin and 
brass column used in our experiments. The surfaces were 
scanned on 800 [im by 800 fim, with 2 /im steps; vertical 
resolution is 100 nm. The horizontal axis on the graphes is 
the symmetry axis of the columns. 



ness 7 /zm with undulations of wavelength 100 /im and 
depth 25 /im. 

In order to compare precisely rheological properties of 
granular assemblies to solid friction properties, we built 
a special device (called "the slider", see Fig. ^ inset) de- 
signed to study the tribology of bead/ wall contact. This 
device is set to apply a constant normal load (Fn = 2 N) 
on three steel beads sliding vertically on the cylinder's 
wall. Then, the dynamical evolution of the resistance 
force encountered by the piston pushing a granular ma- 
terial can be compared to the slider's friction resistance 
driven in the same conditions. 

Two granular systems will be studied. First a model 
assembly of monodisperse steel spheres and second a 
polydisperse assembly of rough glass beads. The largest 
part of this report is devoted to the monodisperse assem- 
bly 



III. MONODISPERSE STEEL BEADS 

We now report on the simplest situation i.e. monodis- 
perse low frictional steel beads in a duralumin or brass 
cylinder. We observe two distinct dynamical stationary 
regimes (Fig. [3J: for high driving velocities, the motion is 
characterized by a steady-sliding and a constant pushing 
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FIG. 3: Resistance force to pushing of a height H = 2.15D 
(380 g) of steel beads in the duralumin cylinder vs. dis- 
placement dstage of the translation stage, a: in the stick-slip 
regime (V = 30 nm.s" 1 ). b: in the steady-sliding regime 
(V = 100 ^m.s- 1 ). 



FIG. 4: Evolution of the resistance force in the steady-sliding 
regime for two different initial packings: a dense one (65%) 
and a loose one (59%); the full line is for the initially dense 
packing, the dotted line is for the initially loose packing, a) 
V = 20 /xm.s -1 ; b) V = 100 /Ltm.s -1 . 



force; for low velocities, the system undergoes a dynamic 
instability characterized by a stick-slip motion. 

We first notice that as far as averaged packing frac- 
tion and pushing force are concerned, whatever the ini- 
tial state of preparation is, the same stationary regime 
is reached (Fig. @J. When a dense packing (D — 
65.0 ± 0.5%) is prepared by rain-filling, the force in the 
steady sliding regime increases up to a maximum, then 
decreases slowly for a piston displacement of about 3 mm. 
At this point where a stationary value F of the force is 
attained, the average packing fraction v is 62.5 ± 0.5%. 
In the case of an initially loose packing [y = 59.0 ±0.5%) 
prepared by using an inner cylinder slowly removed af- 
ter filling, the force in the steady sliding regime increases 
monotonically for a 3 mm displacement before the same 
stationary value F of force is reached. This state is also 
characterized by a 62.5% packing fraction. For all prepa- 
rations, we always get a stationary regime characterized 
by the same pushing force F, which depends on driving 
velocity V", relative humidity \ and packing's height H , 
as well as the same packing fraction [v = 62.5±0.5%), in- 
dependent of V, x an d H. The stick-slip regime displays 
a similar phenomenology: the piling reaches a stationary 
packing fraction of v = 62.5 ± 0.5%; in the transitory 
regime, the pushing force for dense packing displays a 
maximum and for a loose packing shows a monotonous 



increase. In section lill A ll we propose an explanation for 
the phenomenology of the transitory regimes. 

In the following, we present the phenomenology ob- 
served in the steady-sliding and the stick-slip regime, 
when the stationary state is attained. For each regime, 
we first focus on the mechanical properties and propose 
a model accounting for friction at the walls and bulk 
properties; we then study and analyse the rheological 
properties. In a third section, we point out ambiguity of 
dependence on relative humidity. We then analyse our re- 
sults in a standard solid on solid friction framework (the 
Dieterich-Ruina model) . Finally, we study the transition 
from steady-sliding to stick-slip. 

A. Steady-sliding 

For a given height of beads in the column, at given 
driving velocity V and relative humidity x, the force in 
the steady sliding regime is constant at about 2% (Fig. 

and its distribution around mean value F is nearly 
gaussian. 

Next, we study the mean pushing force F behavior as 
a function of the packing height, the driving velocity and 
the relative humidity. 
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FIG. 5: a: Distribution of the resistance force value for every 
fim during 5 mm of sliding in the steady-sliding regime of Fig. 
E|at V — 100 (im.s -1 . b: Distribution of the maximum and 
minimum resistance forces F max and F m i n in the stick-slip 
regime of Fig. [3]at V = 30 nm.s -1 during 5 mm of sliding. 

1. Mechanical properties 

For a vertically pushed granular assembly, the driving 
force exerted by the piston is screened due to friction with 
the walls. To evaluate this effect, the mean resistance 
force F in the steady-sliding regime is measured as a 
function of the packing height (see Fig. |BJ|. 

The resistance force F increases very rapidly with the 
packing's height H. Following the standard Janssen 
screening picture, this strong resistance to motion is due 
to the leaning of the granular material on the walls cre- 
ated by the horizontal redirection of vertical stress in 
association with solid friction at the sidewalls. It means 
that we may relate horizontal and vertical stresses aver- 
aged on a slice at height z by an effective relation: 

a zz {z) = Ka rr (z) (1) 

where K is called Janssen's redirection constant. At the 
walls, we suppose a sliding of the granular material at a 
velocity Vo (the driving velocity); the shearing stress is 
then 

<7 rz (z) = efi d (V Q )a rr (z) (2) 

where HdiVo) is the dynamic coefficient of friction be- 
tween the beads and the cylinder's wall at a velocity Vq. 
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FIG. 6: Mean resistance force in the steady sliding regime 
as a function of the height H of the packing scaled by the 
column diameter D, for steel beads in the duralumin column, 
for V — 16 /im.s -1 (filled squares) and V = 100 fj.rn.s~ 1 (open 
circles). The line and the dotted line are fits by eq. The 
dashed line is the hydrostatic curve. 

The constant e = ±1 is introduced in order to differen- 
tiate between pushing and pulling experiments. When 
e = +1, the granular material is moving upwards and 
friction is fully mobilized downwards (our pushing exper- 
iment) and when e = — 1, the granular is moving down- 
wards and friction is fully mobilized upwards. 

The force F £ exerted by the grains on the piston can 
be derived from equilibrium equations for all slices, thus 
we obtain: 

F e = ggXnR 2 x e(exp(e-) - 1) (3) 

where g is the mass density of the granular material, R is 
the cylinder radius and g the acceleration of gravity. The 
length A = i?/2if^ ( j(Vb) is the effective screening length. 

It is easily seen from (0 that when e = +1, any slight 
change in (i or K is exponentially amplified with a dras- 
tic influence on the pushing force F. In the steady state 
regime, the experimental data obtained for a given push- 
ing velocity V can be fitted by relation © by adjusting 
only one parameter i.e. 

P+1 =Kx fi d (V Q ) (4) 

where Hd(Vo) is the dynamic coefficient of friction at ve- 
locity Vb. 

We see on Fig. [H]that the data are well fitted by eq. 
and, for a relative humidity \ — 45%, for steel beads in 
the duralumin cylinder, we obtain p +1 = 0.140 ± 0.001 
at V up = 16 ^m.s -1 and p + i = 0.146 ± 0.001 at V up = 
100 Atm.s" 1 . 

We are now able to propose an interpretation of the 
features observed for different initial packing fractions 
within Janssen's framework. Since Janssen's coefficient 
K is an increasing function of the packing fraction, as 
it was shown previously |l5l Il7| . the resistance force 
of an initially dense packing displays a maximum, as a 
dense packing leans more efficiently on walls than a loose 
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FIG. 7: Resistance force to pushing (filled squares) and to 
pulling (open squares) in the steady-sliding regime at V = 
16 /im.s -1 as a function of the height H of the packing scaled 
by the column diameter D, for steel beads in the duralumin 
column. The lines are the fit with eq. © of the resistance 
force to pushing and its prediction for the pulling situation; 
the dotted lines are the fit with eq. of the resistance force 
to pulling and its prediction for the pushing situation. 



F + i to F-i. In Fig. 0the pushing force F_i is measured 
for different packing heights H. The fit of experimental 
results with eq. © gives p_i(16 /<m.s -1 ) = 0.156±0.002 
which is 10% larger than p +1 (16 /im.s -1 ). This differ- 
ence, though small, can be observed out of uncertainties, 
and is systematic. It cannot be due to a slight change 
in compacity v as from relation AK/K w 5AD/D, we 
would expect a 2% variation in compacity between the 
pushing and the pulling experiment, which would be ob- 
served; we actually measured Av/v = 0± 1%. According 
to Janssen's picture, this would imply that vertical stress 
redirection is more efficient in the downward pulling sit- 
uation. We believe this is a clear evidence of a granular 
structuration effects but its also shows that this effect is 
not dominant: it affects only 10% of the average mechan- 
ical parameter K. 

We have shown in another report that in a Janssen ex- 
periment, an isotropic homogeneous elastic material can 
also be characterized by stress redirection properties |15| . 
For a granular column at high depth an effective relation: 



one; the pushing force then decreases while the packing 
loosens (i.e. while K decreases). On the other hand, for 
an initially loose packing the resisting force in the transi- 
tory regime will increase while the packing densities (i.e. 
while K increases). 

For an initially dense packing, we extract a parameter 
pmax f rom a f t ue pushing force maximum value 
in the transitory regime with formula ©: we obtain 
p max _ o.l8 at V up = 16 /mr.s -1 and p™? x » 0.18 at 
V up — 100 /im.s -1 ; as this experiment was not repeated 
sufficiently, we do not have uncertainties on these values. 
Nevertheless, we notice that the p™^ x value is roughly 
25% higher than the p + \ stationary value, whatever the 
velocity is. These results are consistent with a K depen- 
dence on compacity v derived in 0] for an assembly of 
monodisperse glass beads, i.e. AK/K w 5Ai>/i>; for ini- 
tial compacity 65% and steady-state compacity 62.5%, 
this empirical formula actually leads to AK/K « 0.2. 
Therefore, in the Janssen framework, we are led to at- 
tribute the difference between maximum force and sta- 
tionary force to a difference in stress redirection (i.e. in 
K) due to a difference in compacity. The force history 
in the transitory regime would then just reflect the com- 
pacity history. 

As a check of consistency, we perform the following 
dynamical experiment. First, the granular column is 
pushed upwards in order to mobilize the friction forces 
downwards and far enough to reach the steady state com- 
pacity. Starting from this situation, the friction forces at 
the walls are reversed by moving the piston downwards at 
a constant velocity Vdown = 16 /im.s -1 , until a station- 
ary regime is attained. Note that this stationary regime is 
characterized by the same compacity v 62.5% as in the 
pushing situation. Following relation 0, this procedure 
would imply a change of e from 1 to — 1 , and consequently, 
the dynamical force on the piston should decrease from 



a zz — K et cr rr (5) 

is obtained with 

K et = v p /{l - i> p ) (6) 

where v p is the material's Poisson coefficient. This actu- 
ally leads to a curve similar to Janssen's saturation curve 
as long as friction at the walls is small (typically less than 
0.5). 

In order to get the isotropic homogeneous elasticity pre- 
diction for the pushing experiment, we perform a series of 
numerical simulations using Finite Element Method • 
The column is modelled as an isotropic elastic medium. 
We vary the friction at the walls /id, the Young mod- 
ulus E and the Poisson coefficient v p . We impose a 
rigid, either perfectly stick or perfectly slip bottom. We 
find no appreciable difference between these two previous 
cases. The condition a rz = —[id^rr is imposed every- 
where at the walls (for the pulling situation, we impose 
<j rz = +/irfcr rr ). The cylinder is modelled as a duralumin 
elastic medium. As long as the Young modulus E of the 
elastic medium is less than 500 MPa, which is usually 
the case for granular media, we find no dependence of 
the results on E. We verified that in all the simulations 
we performed, there is no traction in the elastic medium, 
so that this could be a fair modelling for a granular ma- 
terial. 

We find no appreciable difference between the elastic 
prediction (Fig. |SJ and the curve given by eq. © with 
K = K e i. Therefore, regarding the dependence of the 
stationary state force F on the height of beads, our sys- 
tem cannot be distinguished from an elastic medium. 

Note that finite element simulations show that the 
presence of a rigid bottom implies that the effective 
Janssen's parameter K e ff extracted from Janssen's scal- 
ing for the pulling situation is higher than K e \ |15| . 
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Weight of elastic medium (N) 

FIG. 8: Comparison of the resistance force to pushing sim- 
ulated for a homogeneous isotropic elastic medium (squares) 
of Poisson coefficient v v = 0.45 and Young modulus E = 100 
MPa in a duralumin cylinder, with coefficient of friction 
[i = 0.2 at the walls, to the curve obtained with eq. © 
with Janssen coefficient K = K e i — v v j(\ — v v ) = 0.82 (full 
line). 



whereas for the pushing K e ff « K e i (as can be seen 
on Fig. |SJ the fit of the elastic curve with K = K e i is 
good) . Details on this numerical work will be given else- 
where Actually, if we adjust the elastic predictions 
for pushing and pulling experiments with an clastic ma- 
terial of Poisson coefficient v p = 0.45, eq. © yields a 
Janssen's constant K e ff for the pushing which is about 
3% lower than K e f f for the pulling. This is qualitatively 
(though not quantitatively) in agreement with the exper- 
imental results. Then isotropic elasticity can be a good 
framework only if we neglect the existence of bulk struc- 
turation effects inducing differences in the effective Pois- 
son coefficient of the material between the pulling and the 
pushing. Note that in this case, an isotropic modelling 
of the granular material is somehow questionable. 

In a previous study 01 > ^ was found that the Janssen 
picture has a general tendency to slightly underestimate 
the stress below a granular column for a homogeneous 
packing of glass beads. We have showed in another re- 
port that it is no more true as the friction at walls is 
very well controlled. Therefore, this model, though ele- 
mentary, seems a fair base for analysis and provides an 
analytical expression from which constitutive rheological 
parameters can be extracted. A central question is still 
that the fitting parameters p = K x fid extracted from 
the model does not allow to distinguish between fid and 
K separately. 

In the following, we will show experiments which aim 
at sorting out the relative contributions of wall-bead in- 
teractions (i.e. fid) and bulk properties (i.e. K) which 
have an influence on the rheological properties when ve- 
locity and relative humidity are changed. 
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FIG. 9: a) Mean force F in the steady sliding regime for height 
H — 2.15D (380 g) of steel beads in the duralumin cylinder, 
for various relative humidities X- 72% (diamonds), 66% (open 
inverted triangles), 53% (triangles), 40% (open circles), and 
< 3% (squares) . b) Coefficient Kfid extracted from a fit of F 
with formula J3J. 



2. Rheological properties 

The mean resistance force F in the steady-sliding 
regime increases strongly when the velocity and the rel- 
ative humidity are increased. In the duralumin col- 
umn (Fig. EJ, the resistance force F at velocity V — 
100 ^m.s -1 is 60% higher for relative humidity \ = 73% 
than for \ = 40%; near the transition to the stick-slip 
regime (V — 1 /im.s -1 ), dependence of F on \ is much 
less important, and there is actually no difference in F 
values for \ varying between 40% and 73%. Another im- 
portant feature is that the increase of F with velocity is 
sharper when \ is higher. In the brass column (Fig. I10[) . 
the mean resistance force F also increases with relative 
humidity \i but now the variation coefficient with veloc- 
ity (the slope on Fig. I1U|) does not seem to depend on \- 
Moreover, for a velocity V — 100 /zm.s , the F increase 
is only 15% from dry to humid (\ — 90%) atmosphere. 

The fit of F with Janssen's formula adapted to the 
pushing case of eq. © gives parameter Kfid(V) for dif- 
ferent relative humidities \- The fundamental differences 
of behavior in the steady sliding regime between the du- 
ralumin column and the brass column for a same granular 
material, at same density, suggest that the variations in 
X and velocity have an influence mainly on the coefficient 
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<3% 


0.0065 ± 0.0006 


40% 


0.004 ± 0.0005 


53% 


0.013 ±0.001 


66% 


0.021 ± 0.0005 


72% 


0.028 ±0.0015 



TABLE I: Parameter K(a—b) extracted from a fit of Kfid with 
formula fid(V) = fio + {a — b) ln( V/Vq) for different relative 
humidities \ m the duralumin column. 
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FIG. 10: a) Mean force F in the steady sliding regime for 
height H = 2.15D (380 g) of steel beads in the brass cylinder, 
for various relative humidities \- 90% (squares), 53% (open 
circles), and < 3% (triangles), b) Coefficient Kfid extracted 
from a fit of F with formula . 



FIG. 11: Dynamical coefficient of friction fid extracted from 
a fit of F with Janssen's formula adapted to pushing (open 
symbols) for a height H — 2.15D (380 g) of steel beads in the 
duralumin cylinder, and fid measured in a solid on solid ex- 
periment (filled symbols), for two relative humidities \- 40% 
(circles) and < 3% (squares). 



of friction fid at the walls, and little on mechanical prop- 
erties of the granular material (i.e. K). Therefore, in the 
following analysis, we will consider if as a constant at 
first order of approximation. 



Steel/ duralumin 

Fig. |5|shows that dynamical parameter K fid is globally 
less important in a dry atmosphere than in the ambient 
atmosphere. In ambient atmosphere (40% < x < 73%), 
a change in \ seems to change only the variation coef- 
ficient of Kfid with velocity: i.e. the increase of K fid 
with velocity is sharper when \ is higher. This phe- 
nomenon could be interpreted as a viscous contribution 
of the water condensed at the contacts. Note that this 
phenomenology contrasts with the observations of Riedo 
et al. || who find, for solid on solid nanoscopic slid- 
ing friction measurements, that in all the systems they 
study, of various wettability, the effect of a humidity in- 
crease is to add a negative value to the coefficient of fric- 
tion force dependence on logarithm of velocity. Table |I] 
shows the results of a rough logarithmic fit of data with 
fJ>d(V) = Ho ± (a — 6)ln(V/Vo). We use the parameter 
(a — b) as a standard reference to the Dieterich-Ruina 
model [Ulli. 



We observe that the increase rate K(a — b) of K fid 
with the logarithm of velocity is multiplied by 7 when x 
increases from 40% to 72%. 

In order to compare these results with those obtained 
in a solid on solid friction experiment, we need to evalu- 
ate Janssen's parameter K. We extract K from a clas- 
sical Janssen experiment |l5j and we obtain: Kfi s = 
0.184 ±0.002 for a mean aging time of 40 seconds. From 
the measurement of static friction, using the sliding angle 
of a three bead tripod, we extract K = 1.02 ± 0.07. 
We suppose K is constant and is not affected by any 
variation in relative humidity nor velocity. The dynami- 
cal coefficient of friction obtained in a solid on solid fric- 
tion experiment with the slider (Fig. ^ inset) is then 
compared on Fig. ^]with the one extracted from a fit 
of the resistance force with formula 10 • Note that the 
uncertainty on fid extracted from formula @ is mainly 
systematic (coming from uncertainties on K, which come 
themselves from uncertainty on static friction measure- 
ments). We observe that the dynamical coefficients of 
friction measured in a solid on solid experiment increase 
slightly more strongly with velocity than the one ex- 
tracted from the resistance force of the granular mate- 
rial. But importantly, the increase of fid with relative 
humidity x. is the same in both cases. 
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relative humidity \ 


K(a - b) 


<3% 
53% 
90% 


0.007 ±0.001 
0.007 ±0.001 
0.0085 ± 0.0015 



TABLE II: Parameter K{a — b) extracted from a fit of Kfid 
with formula fid(V) = fj,o + {a.—b) ln(V/Vb) for different relative 
humidities x m the brass column. 



In this analysis we assumed a constant value for Janssen's 
coefficient K. We can also imagine a slight decrease of 
K with increasing velocity so that the solid on solid fj,^ 
would match the one extracted from the pushing exper- 
iment. However, we will remain prudent as the mean 
pressure contact on beads in the solid on solid experi- 
ment is much higher than in the granular column, which 
may cause quantitative differences between friction prop- 
erties. More precisely, at the bottom of the column, for 
a pushing force F = 10 N, the mean contact force per 
bead is about 20 mN; a Hertz contact would give mean 
contact pressure p ~ 300 MPa. For a contact force per 
bead of 2 N (the slider case) , a Hertz contact would give 
p ~ 1300 MPa. 



Steel/brass 

In the brass column, we observe the same phenomenol- 
ogy as in the duralumin one, i.e. an increase of K\id with 
velocity and relative humidity x- The effect of a change 
in x is yet less important than in the duralumin column, 
as the coefficient Kfid increases only by 15% at veloc- 
ity V = 100 /im.s -1 when \ increases from 0% to 90%. 
A major difference with the duralumin column is that 
in the brass column, a variation in x seems to induce a 
variation of friction but hardly affects the coefficient of 
variation with velocity (the slope). Table ITU shows the 
results of a logarithmic fit of the data with the function: 
HdiV) = Mo ± (a — b)ln(V/Vo)- We observe indeed that 
the slope K(a — 6) of K[id with the logarithm of velocity 
does not practically vary with relative humidity x, and 
is much weaker than in the duralumin column. 

The comparison with the slider, and the differences 
between the duralumin and bras columns, suggest that 
the rheological properties in the steady-sliding regime are 
dominated by the friction properties at the walls. We now 
investigate the stick-slip regime. 



B. Stick-slip regime 

In the stick-slip regime, for a given height of beads in 
the column and at given driving velocity V and relative 
humidity Xi the maximum force before a slip and the 
minimum force at the end of a slip are constant within 
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FIG. 12: a: Force vs. time for a slip event in the stick-slip 
regime at V = 1 fim.s -1 for a height H = 2.15D (380 g) of 
steel beads in the duralumin column, b: Force vs. velocity of 
the granular material for the slip event of Fig. a 



2% (Fig. |SJd). Their distribution around mean values 
F max and F mi „ is nearly gaussian. 

The dynamical evolution F(t) of the pushing force dur- 
ing the slip phase can be translated into a function of 
the granular material instantaneous velocity. Motion of 
the granular material at velocity V induces a variation 
dF = fc(Vo — V)dt of resistance force during time dt, 
where k is the force probe stiffness and Vo the driving 
velocity; we thus obtain V(t) = Vq — dF/kdt. We see 
(Fig. I12|) that the acceleration and deceleration phases 
are nearly symmetric. This phenomenology is currently 
observed in solid friction experiments. This is in contrast 
with previous plane shearing experiments of granular ma- 
terials @ in which the deceleration phase would occur in 
two steps, first fast and then slow. However, a major dif- 
ference is that in the latter case, the displacements during 
slippage are millimetric (a grain size) whereas in our ex- 
periment, they are micrometric (30 /mi on Fig. I12f> . As 
a consequence, the evolution of forces observed during a 
slip seems to reflect the slippage of grain contacts at the 
walls. The maximum velocities obtained during the slip 
phase are of a few hundred /zm.s . We did not study sys- 
tematically the variation of instantaneous velocity with 
bead height nor with driving velocity. 

In the following, we report on the F max and F min prop- 
erties as a function of the packing height, the driving 
velocity and the relative humidity. 
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FIG. 13: Mean maximum and minimum resistance forces 
Fmax (filled triangles) and Fmin (open circles) in the stick- 
slip regime as a function of the height H of the packing scaled 
by the column diameter D, for steel beads in the duralumin 
column, for stick-time t 3 u c k = 40 s. The full and dotted lines 
are fit by eq. JSJ. The dashed line is the hydrostatic curve. 

1. Mechanical properties 

In the stick-slip domain, the mean maximum and min- 
imum resistance forces F max and F m i n are measured as 
functions of the packing height (see Fig. 1131) , but now we 
choose to perform our experiments not at a constant ve- 
locity V, but at constant stick time tstick and we explain 
why in the following. 

Experimentally we find that the resistance forces 
F m ax and F m i n increase exponentially with the packing's 
height H. The data can actually be fitted by a formula 
of type eq. J3J), but the parameter p + i extracted may be 
different and should not depend on /id(Vb) as the granu- 
lar material is not sliding during stick, and does not slide 
at driving velocity Vo during slip. 

The model can actually be simply modified for this sit- 
uation. In the stick-slip regime, when the resistance force 
F(t) is equal to F max , the material just starts slipping, 
which means that static friction is fully mobilized at this 
instant for each contact. So we can write a rz — [i s cr rr ev- 
erywhere at the walls, where [i s is static friction between 
beads and the cylinder's walls. Static friction coefficients 
are known to evolve in time |2l| . we then have to in- 
clude aging effects. As the granular material was at rest 
in the column since the previous slip, aging time in this 
situation is the time of stick: 

, _ F max F m z n . . 

"stick TTT V ' / 

kVo 

So the final formula for F max in the Janssen framework 
is: 

TT 

F e = ggXnR 2 x e(exp(e— ) - 1) (8) 
A 

where the screening length is now A = R/2K [i, s {t stick)- 
Here again, e = +1 for a pushing experiment (i.e. 
F m ax = F +1 ), e = — 1 for a Janssen experiment. The 



FIG. 14: Sketch of the method for measuring forces at the 
onset of slipping in the pushing and pulling situations, (a) 
refers to the maximum force before a slip in the stationary 
stick-slip regime in the pushing experiment; (b) refers to the 
minimum force before the first slip in the pulling experiment 
after the material have been pushed over 5 mm; (c) refers to 
the minimum force before a slip in the stationary stick-slip 
regime in the pulling experiment. 

choice of imposing a constant stick-time t st i C k instead of 
constant driving velocity Vq for the experiments in the 
stick-slip regime is now justified by the dependence of /x s 
on t S uck, whereas [id depends on Vq. 

In the steady state regime, the experimental data ob- 
tained for a given stick-time t st i C k can be fitted by re- 
lation (jSJ) by adjusting only one parameter i.e. p+i = 
K x [i s (t stick) where fj, s (t stick) is the static coefficient of 
friction for an aging time t stick- For the data of Fig. 
IT31 obtained for a relative humidity \ — 40% with a 
mean stick time of 40 s, we find p+i = 0.184 ± 0.002 
from F max (F +1 in eq. JSJ. Note that, though friction 
is not fully mobilized when F = the fit of F m i n 

with eq. © is good (Fig. ITBjl : we obtain a parame- 
ter p + i = 0.148 ± 0.002 which should not be related to 
Vs(t stick) in that case. 

As a check of consistency, we perform the following 
dynamical experiment in the stick-slip regime. First, the 
granular column is pushed upwards in order to mobilize 
the friction forces downwards and to reach the steady 
state compacity; F ma x (referred to as (a) on Fig. I14|) 
is then measured for constant stick time t s tick- Starting 
from this situation, the friction forces are reversed at 
the walls by moving the piston downwards. We then 
measure the force F m in of the first slip (referred to as 
(b) on Fig. I14|) after the same time t s tick- Note that 
the roles of F m i n and F max are inverted in the pulling 
experiment: F m in is the minimum force before a slip, and 
Fmax is the maximum force at the end of a slip. We also 
measure mean minimum force F m i n (referred to as (c) 
on Fig. I14f) when a steady stick-slip regime is reached in 
the pulling experiment. Note that this stationary regime 
is characterized by the same compacity v w 62.5% as in 
the pushing situation. 

Following relation (jHJ, this procedure should imply a 
change of e from 1 to — 1 , and consequently, the force at 
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Filling mass (g) 

FIG. 15: a: Mean maximum resistance force measured in the 
stick-slip regime as a function of the height H of the packing 
scaled by the column diameter D, for steel beads in the du- 
ralumin column, for stick-time t stic k = 40 s; the line is a fit 
with eq. JHJ. b: Measured mass vs. filling mass for a classical 
Janssen experiment obtained in two different ways: we plot 
the mass of the first slipping event (squares), and the mean 
minimum mass in the stationary stick-slip for the pulling situ- 
ation (triangles); the line is the Janssen curve predicted with 
the parameter extracted from the pushing experiment of Fig. 
ITSl . 



p+i(40 s). We believe this is a clear evidence of a gran- 
ular structuration effects. We interpret the differences 
in p_i values by saying that the packing has first been 
structured to resist the pushing; when we start pulling, 
structure is not efficient to resist the pulling, and it re- 
sults in a lower Janssen's constant (therefore lower P-i) 
for the first slip event; when we continue pulling, the 
packing gets structured to resist pulling and leans more 
efficiently on the walls, i.e. Janssen's constant (and p-i) 
increases. Moreover, steady-state structuration is more 
efficient for the pulling than for the pushing, which means 
that both structurations may be different and could re- 
flect the symmetry breaking due to gravity. 

The results of the simulations of an elastic medium in 
a column presented in the precedent section for a steady- 
sliding also carry on to the stick-slip situation. If we im- 
pose numerically everywhere at the walls, the relations 
a rz = (J> s cr rr or <j rz = fi d a rr , it would merely correspond 
to a change in the name of the friction coefficient. The 
equilibrium equations are the same for the static and 
steady dynamic cases. The important common point is 
that friction forces are fully mobilized everywhere at the 
walls. Therefore, as for the steady-sliding, the data can 
be fitted by the elastic predictions, and the 3% differ- 
ence in Janssen's parameter between the pulling and the 
pushing of an elastic medium of constant Poisson ratio is 
again qualitatively (though not quantitatively) in agree- 
ment with the experimental results. 

In the following, as we study the steady-state regime, 
we will try to distinguish between wall-bead interactions 
(i.e. fj, s ) and the bulk properties (i.e. K) as they influ- 
ence the rheological properties when velocity and relative 
humidity are changed. 



2. Rheological properties 



the onset of slipping on the piston should decrease from 
F +1 to F^ 1 . As we observed no noticeable variation of 
compacity between the pushing and the pulling, we ex- 
pect Janssen's parameter K value to be unchanged. In 
Fig. E| the pulling force F _1 is measured for different 
packing heights H, for the first slipping event and in the 
steady stick-slip regime. The fit of experimental results 
with eq. © gives p_i(40 s) = 0.230 ± 0.004 for sta- 
tionary stick-slip which is 20% larger than p + i(40 s). It 
cannot be due to a slight change in average compacity v 
as from relation AK/K « 5Av/v, we would expect a 4% 
variation in compacity, which would have been observed; 
let us recall that we actually measured a compacity vari- 
ation between the pushing and the pulling experiment 
Av/v = 0± 1%. 

As in the steady-sliding experiment, according to 
Janssen's picture, this would then imply that vertical 
stress redirection is more efficient in the pulling situa- 
tion; furthermore, this effect is enhanced in the stick-slip 
regime. For the first slip event in the pulling experiment, 
we find p_i(40 s) = 0.16 ±0.015 which is now lower than 



The mean maximum force F max increases strongly 
when the driving velocity Vo is decreased and the rel- 
ative humidity x is increased (Fig. ITTik. IT7b. IT%b). Fur- 
thermore, F max apparent divergence when Vq decreases 
is enhanced when \ increases. The mean minimum force 
Fmin is almost independent of the driving velocity, and 
increases when the relative humidity is increased (Fig. 
116b. 117b. 118b). For each stick-slip event, an aging time 
(the stick time) is measured: t stlck = {F max -F mm ) / (kV) . 
This time is associated to a parameter K ^ s it stick) ex- 
tracted from the fit of the maximum force F max of this 
stick-slip event with equation (jSJ. In the following, the 
results of this analysis for duralumin and brass columns 
are compared and discussed. 



Steel/ duralumin 

For steel beads in the duralumin column, the param- 
eter K[i s (Fig. 116b . 117b ) increases roughly logarithmi- 
cally with aging time, and increases strongly with rel- 
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FIG. 16: a: Minimum and maximum resistance forces F m in 
(circles) and F max (squares) vs. driving velocity in the stick- 
slip regime for a height H = 1.4D (250 g) of steel beads 
in the duralumin cylinder, for relative humidities x = 90% 
(filled symbols) and x = 45% (open symbols), b: Kfi s as a 
function of stick time t st i C k for a height H — IAD (250 g) of 
steel beads in the duralumin cylinder, for relative humidities 
X = 90% (filled squares) and x = 45% (open circles). 



FIG. 17: a: Minimum and maximum resistance forces F m in 
(circles) and F max (squares) vs. driving velocity in the stick- 
slip regime for a height H — 2.15D (380 g) of steel beads 
in the duralumin cylinder, for relative humidities x = 45% 
(open symbols) and x < 3% (filled symbols), b: K[i s as a 
function of stick time t stic k for a height H — 2.15D (380 g) of 
steel beads in the duralumin cylinder, for relative humidities 
X = 45% (open circles) and x < 3% (filled triangles) . 



ative humidity x- Its value increases by 50% when x 
increases from 3% to 45%, and doubles when x increases 
from 45% to 90%. The parameter Kfi s does not depend 
only globally on x, but its logarithmic aging rate Kb 
with stick time t st i C k depends strongly on x- The param- 
eters Kb extracted from a logarithmic fit of K[i s (t) with 
/i s (t) = /io + &ln(i/ro) are given in table Hill for different 
relative humidities. 

Note that magnitude of Kb, as K is of the order 1, is con- 
sistent with many previous reports on solid friction 0,0. 
An important feature is that the logarithmic aging rate 
appears to be 4 times higher in humid atmosphere than 
in dry atmosphere. 

The logarithmic evolution is not perfect: for aging time 
larger than 3000 s, we observe that K/i s increases more 
sharply with time; for H = 2A5D of steel beads and 
X = 45%, Kb varies from 0.009 for times t s u c k < 3000 s 
to about 0.025 in the last time decade (Fig. 117b). This 
phenomenon was observed for all relative humidities x- 
It is similar to what was reported previously by Losert 
et al. 0] for plane shearing of glass beads. 
It is tempting to link this enhanced aging behavior to 
the granular material slow restructuring effects since one 



relative humidity x 


Kb 


<3% 
45% 
90% 


0.006 ± 0.001 
0.013 ± 0.003 
0.023 ± 0.001 



TABLE III: Parameter Kb extracted from a fit of Kfi„(t) with 
fis(t) = /io + b\n(t/ro) for different relative humidities Xi f° r 
a height H = IAD (250 g) and a height H = 2.15D (380 g) 
of steel beads in the duralumin column. 



may think of an aging of the granular structure leading 
to a better stress redirection at the walls. 

However, at this point, it is not possible to differentiate 
between wall and structure properties. We tried to mea- 
sure independently /i s aging properties in a solid on solid 
experiment, but we did not obtain satisfactory statistics 
to conclude. Therefore, we need to compare the results 
obtained in the duralumin column to the one obtained in 
the brass column. 
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FIG. 18: a: Minimum and maximum resistance forces Fmm 
(circles) and F max (squares) vs. driving velocity in the stick- 
slip regime for a height H = 1.7 D (300 g) of steel beads 
in the brass cylinder, for relative humidities x = 90% (open 
symbols) and x < 3% (filled symbols), b: K[i s as a function 
of stick time t st i C k for a height H — 1.7D (300 g) of steel 
beads in the brass cylinder, for relative humidities \ = 90% 
(filled squares) and \ < 3% (filled triangles), and for a height 
H = 2.15D (380 g) for relative humidities x = 64% (open 
circles) and x < 3% (open triangles). 



Steel/brass 

For steel beads in the brass column, the parame- 
ter K fj, s extracted from the experimental measurements 
(Fig. 118b ) also increases logarithmically with aging time. 
However, the increase with relative humidity \ is much 
less important than in the duralumin column since K/i s 
increase is only 50% when increasing \ from dry atmo- 
sphere to 90% humidity, whereas the increase is more 
than 150% in the duralumin column for the same vari- 
ation of x- Moreover, the K[i s increase with stick time 
seems now perfectly logarithmic over 4 decades of time 
variation. 

As in the duralumin column, K[l s increase rate with 
stick time t st i C k depends on \. The parameters Kb ex- 
tracted from a logarithmic fit of Kfj, s (t) with fJ, s (t) — 

+ &ln(t/ro) are given in table II VI Logarithmic aging 
rate is now 2.5 times greater in humid atmosphere than 
in dry atmosphere (it was 4 times higher in the duralumin 
column) . 



relative humidity x 


Kb 


<3% 
64% 
90% 


0.006 ± 0.001 
0.011 ± 0.0005 
0.014 ± 0.0005 



TABLE IV: Parameter Kb extracted from a fit of Kfi s (t) with 
jUs(t) = Mo + bln(t/ro) for different relative humidities x, f° r 
a height H = 1.7D (300 g) and a height H = 2.15D (380 g) 
of steel beads in the brass column. 



Comments 

The only differences between both systems (duralumin 
and brass column) being the contacts at the walls, as 
for the steady-sliding regime, the differences observed in 
the rheological behavior suggest that the aging of K[i s 
is mainly an aging of friction coefficient fj, s at the walls 
and that humidity affects principally the contacts at the 
walls. 

Relative humidity x seems to have two effects in the 
duralumin column: first, the aging rate is higher when 
X is higher, second, the friction level for short times is 
higher when x is higher. 

In the brass column there seems to be an effect of x only 
on the aging rate: the Kfj, s (t 8 ti c k) curves for different x 
cross for t — 0.1 ms. The logarithmic aging rate coeffi- 
cient gets comparable values in both columns, but it is 
higher in the duralumin column (Kb — 0.023) in humid 
atmosphere than in the brass column (Kb = 0.014). 

Our observations are consistent with recent aging ex- 
periments in granular media |22j,|23j and solid on solid ag- 
ing experiments 4] , who found logarithmic aging of static 
friction enhanced by an increase in relative humidity x- 
These results have been interpreted by the dynamics of 
capillary condensation at the contacts [22L |24| , which is 
a thermally activated process: as time goes on, there are 
more and more capillary bridges at the contacts, which 
are responsible for an adhesion force increasing with ag- 
ing time; moreover, condensation goes faster when x is 
higher. In these models, there is no aging in dry at- 
mosphere, as no condensation can occur; this prediction 
agrees with most observations Q. As in our experiment 
the logarithmic aging rate coefficient is not zero in dry 
atmosphere (even though the air flux we impose is prob- 
ably not perfectly dry), we can conclude that capillary 
condensation is not the only source of aging. Another 
source of aging can be creep of contacts Q ■ 

We observed that aging is perfectly logarithmic in the 
brass column. Therefore, we are led to analyze again the 
aging in the duralumin column since the interpretation 
based on aging of internal friction we proposed for the du- 
ralumin column, should also apply to the brass column, 
which is in contradiction with the observed behavior. 

If we suppose that K properties depend only on the 
granular material properties, the acceleration of logarith- 
mic aging in the duralumin column may be not linked to 



aging of bead/bead contacts nor granular slow restruc- 
turations. We may then start from the following exper- 
imental observations: in solid on solid friction experi- 
ments, static friction coefficients /j. s are found to depend 
on the applied shear exerted for a given waiting time 
pi l2rj| , For a higher shear during a given waiting time, 
coefficient fi s is higher, and the logarithmic aging rate 
is higher too 0. In our experiments, during a "stick" 
event, the pushing force increases linearly with time, i.e. 
shear at the walls (and contact pressure) increases. As 
a consequence, aging occurs with a non-constant applied 
shear and the mean shear is more important when the 
aging time is higher. Consequently for high waiting time 
or stick time, i.e. for high mean shear at the walls, the 
aging of the static friction coefficient may be accelerated. 
This phenomenology may be a priori different for differ- 
ent surfaces, which would explain why aging seems per- 
fectly logarithmic in the brass column, and not in the 
duralumin column. 



3. Analysis of Fmax evolution 

Now, the high enhancement of the blocking resistance 
force F m ax with decreasing driving velocity can be sim- 
ply understood within the Janssen framework . We start 
from the experimental observation that the minimum 
force F m i n after a slip does not depend on velocity (at 
constant relative humidity y): Fmin — Fq. Therefore, 
the aging time t st i C k before a slip reads: 

tsuck = (F max -F )/(kV) (9) 

We suppose, as well verified experimentally, that the co- 
efficient of friction fi s at the walls follows a logarithmic 
aging law: 

M s (*) = Mo+&ln(t/r) (10) 

This value of /U s (t) for aging time t sti( ± can be injected 
in formula (JHJ. We then obtain: 

F max cx V- 1 ^ (11) 
with 7 = R/(2KbH) - 1. 

The apparent divergence of F max with decreasing veloc- 
ity, as observed experimentally, is thus more important 
when the height of beads and the logarithmic aging rate 
b of friction at the walls (i.e. relative humidity as b in- 
creases with y) are higher. The observed strong depen- 
dence of F max with V then corresponds to an exponential 
amplification of friction aging at the walls. Elements of 
interpretation along these lines were already given in rcf 

m 

C. Influence of relative humidity: an open problem 

Nevertheless from the whole series of experiments we 
performed on this system we are led to conclude that 
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FIG. 19: a: Maximum and minimum resistance forces to 
pushing Fm ax (stars) and F m i n (squares) of a height H — 
2.15-D (380 g) of steel beads in the duralumin cylinder vs. 
time in the stick-slip regime (V = 50 nm.s -1 ). b: Evolution 
of relative humidity y during the experiment of Fig. a. The 
dotted line separates two phases: 1) forces and x evolutions 
are not correlated; 2) forces and x evolutions are correlated. 
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FIG. 20: Maximum and minimum forces F max (stars) and 
F m in (squares) vs. relative humidity x f° r the experiment of 
Fig. E0 phases 1) and 2) of Fig. EH arc indicated on the 
graph. 



the dependence of the phenomenology with the relative 
humidity x remains sometimes unclear. 

As an illustration of this statement, we report the fol- 
lowing result (Fig. 1191 120fl . In an experiment performed 
in the stick-slip regime and in the ambient atmosphere, 
for a driving velocity V = 50 nm.s -1 , the forces F m j„ and 
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Fmax were found to be independent of x, which was vary- 
ing between 52% and 62%, during the first 30 hours of the 
experiment. We then observed that the F m i n and F max 
variations became suddenly correlated to the \ variations 
and would stay correlated for the next 80 hours. For 
the first 30 hours, x increased from 52% to 62%, while 
F mm and F max were nearly constant: F rnax = 11.6 N, 
F mm = 8.2 N. Then x varied from 62% to 68% in one 
hour: F m i n and F max started to be correlated to x from 
this " triggering event" and for the next 80 hours as if the 
system had suddenly reached some "reversible branch" . 

Such a phenomenology was never obtained when hu- 
midities x = 90% and x < 3% were imposed. When we 
imposed dry atmosphere, the force would reach its sta- 
tionary value in a few minutes, and we never noticed any 
variation in this value for 5 days; when we imposed humid 
atmosphere, the force would reach its stationary value in 
a maximum of 3 hours and we noticed no variation in 
this value for 5 days. 

This may however incite to prudence for results ob- 
tained in ambient atmosphere, for experiments per- 
formed on the granular column as much as for solid on 
solid friction experiments. As an example, we see on Fig. 
EDI that for x = 60%, F max = 11.5 ± 0.2 N when forces 
and x are n °t correlated, whereas F max = 10.3 ± 0.2 N 
when they are correlated. 

Furthermore, this opens the question of the metasta- 
bility of capillary condensation. If variation of F max with 
X is interpreted as variation of fi s due to thermally acti- 
vated capillary condensation at the contacts [22, |24| , it 
seems possible that capillary condensation did not occur 
for hours. If we notice that F max starts to be correlated 
with x when F max constant value is equal to F max on 
the x dependent part of the curve (i.e. F max — 11.8 N 
and x = 68%), another possibility is that there were as 
much capillary bridges condensed at the contacts as for 
X = 68% in the first part of the experiment, i.e. that 
evaporation does not always happen. 



D. Dieterich-Ruina analysis 

A standard phenomenological model accounting for 
static and dynamic properties of solid friction is the so- 
called Dieterich-Ruina model [l^. l20| . In this picture, 
the coefficient of friction is 



This model accounts for the logarithmic aging of static 
friction, and the logarithmic velocity strengthening or 
weakening in the steady-sliding regime with: 



fi = fi + a\n(—) + bin —- 
Vo L) 



where 9 obeys 2t 



<W_ 6V 
dt~ D 



(12) 



(13) 



fio, fl j b, and D are constants dependent of the mate- 
rials in contact. Parameter D is usually interpreted as 
a characteristic length for renewal of contacts and is of 
order of the creep length before sliding |23|. In a sta- 
tionary steady-sliding regime 9 = D/V is interpreted as 
a characteristic time for renewal of contacts. 



d/i s /dln(i) = b 
d^ d /dln(V) = a - b 



(14) 
(15) 



It also accounts for a third usual phenomenon observed 
in the transitory regime when the velocity is changed 
suddenly. In this case there is an abrupt change in fric- 
tion coefficient followed by exponential relaxation to a 
new stationary value. If we change velocity V\ to ve- 
locity V2, the instantaneous change in friction coefficient 
from stationary value Hd(Vi) = /io + (a — b)\a{V\/Vo) 
is +aln(V2/Vi); friction change from this value to new 
stationary value fidi}^) is then — &ln(V2/Vi) during typ- 
ical time t = D/V2, i.e. we get ^(^2) = ^0 + ( a ~ 
b)\n{V 2 /V Q ). 

Now we test the Dieterich-Ruina model with the 
steel/duralumin and steel/brass friction coefficient ex- 
tracted from our data. As we observed logarithmic aging 
and rough logarithmic velocity strengthening, we get pa- 
rameters a and b. We are therefore able to test this model 
predictions with the observation of response to an abrupt 
variation of velocity. 

Experimentally, we never observed any transitory 
state: the force in the steady sliding regime was always 
changed from stationary value F(V\) to stationary value 
F(V2) when changing velocity from V\ to V%. From the 
measured value of Kb = 0.010 for steel/duralumin fric- 
tion and Kb = 0.011 for steel/brass in ambient atmo- 
sphere, we would however expect the force at the begin- 
ning of the transitory regime to be 0.5 N larger (steel 
beads height H = 2.15D) than the stationary force at 
the new velocity. As the natural fluctuations of force are 
around 0.1 N, we would therefore expect to observe this 
transitory regime in F due to change in friction at the 
walls. Now we give an interpretation for why, within the 
Dictrich-Ruina picture, this transitory regime was not 
observed, which yields an upper bound for the material 
contact renewal length Do- 

When velocity is changed from 10 /im.s -1 to 
100 ^m.s- 1 for H = 2.15D (380 g) of steel beads in a 
duralumin cylinder, the change in F stationary value is 
of 1 N. But, due to the finite stiffness of our set-up, this 
change cannot be instantaneous: an increase AF of force 
at constant velocity takes minimal time t — AF/(kV). 
So an increase of 1 N of force at 100 ^tm.s -1 takes time 
t = 0.25 s. This time must be compared to time length 
t of transitory regime: r = Dq/V. Therefore, in the 
framework of the Dieterich-Ruina model, we need r to 
be less than the time of force variation, which means 
Dq < 25 fim . 

Let us now consider a sudden decrease in velocity from 
100 /im.s -1 to 1 fim.s -1 . Now the decre ase of force is 
governed by inertial time: U n = n^m/k which can be 
evaluated around tin = 200 ms. This leads to Do < 400 
nm (so that r < t in ) consistently with the Dieterich- 
Ruina model if no transitory state is observed. 
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FIG. 21: Mean maximum and minimum forces F max (squares, 
dotted line) and F m i n (circles, dashed line) in the stick-slip 
regime, and mean resistance force F (triangles, full line) in 
the steady-sliding regime for a height H = 2.15D (380 g) of 
steel beads in the duralumin cylinder, for relative humidity 
X = 40%. 



We have no precise measure for Dq but this upper limit 
on the value is coherent with the usual interpretation in 
terms of length for contact renewal, as the duralumin 
column mean roughness is 400 nm. Thus, we can neither 
validate nor rule out the Dieterich-Ruina model for our 
system. 



E. Transition mechanism 

For experiments performed at given height of beads 
and relative humidity Xi we observe that the transition 
from stick-slip to steady-sliding occurs for a finite stick- 
slip amplitude (Fig. 121(1 . This suggests a subcritical 
transition but such a transition would be hysteretic. In 
order to verify this feature, we make the following exper- 
iment: we drive the system at constant velocity in the 
steady-sliding regime near the transition until it reaches 
a stationary state, we then decrease continuously the ve- 
locity, without ever stopping the movement, to a veloc- 
ity for which stick-slip used to occur. Note that the re- 
sults reported in the previous sections and on Fig. [2"71 
were obtained for independent experiments on systems 
driven at a unique velocity. We observe that if we impose 
V = 700 nm.s -1 directly, we obtain a stick-slip motion, 
whereas when we decrease velocity from 1.4 (im.s^ 1 (in 
the steady-sliding regime) to V — 700 nm.s -1 the motion 
is steady-sliding. This proves the hysteretic character of 
the transition, which was actually hard to observe as the 
system seems very sensitive to noise. 

This transition is similar to the one observed by Heslot 
et al. |27j in solid on solid experiments; it correspond to 
what they call the " inertial regime" . It is a characteristic 
of systems for which the dynamical coefficient of friction 
increases with velocity. This transition scenario was ex- 
plained by Brockley et al. However they find that 
the force F in the steady-sliding regime at the transition 
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FIG. 22: a: Minimum and maximum resistance forces F m in 
(open circles) and F max (full squares) vs. driving velocity 
in the stick-slip regime for a height H — 2.8D (132 g) of a 
miscellany of glass beads in the PMMA cylinder at relative 
humidity \ = 35%. b: K(X S as a function of stick time t a u c k 
extracted from F max with equation JHJ. 



is of order of F m j„ . We actually find (see Fig. I21|) that F 
is less than F m i n for steel beads in the duralumin cylin- 
der, a feature we do not understand, whereas F rs F m i„ 
for in the brass cylinder. 

Again, the transition we observed is consistent with 
a picture of solid friction sliding instability at the wall. 
This is also consistent with the model we have developed 
in the previous section. 



IV. POLYDISPERSE GLASS BEADS 



After the study of the simplest case, with monodisperse 
low friction beads, we now want to study the effect of 
disorder (friction, polydispersity) on rheology. Therefore, 
we choose to study a miscellany of glass beads using a 
three diameters mixture (1.5 mm, 2 mm and 3 mm) with 
equal volume of each kind in an abraded PMMA cylinder. 
In the following, we present the features obtained in the 
stick-slip regime for this system. 
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FIG. 23: a: Resistance force to pushing of a height H — 2.7 D 
(126 g) of a miscellany of glass beads in the PMMA cylinder 
vs. displacement dstage of the translation stage in the stick- 
slip regime (V = 1500 nm.s -1 ). b: zoom of Fig. 123b . 



FIG. 24: Distribution of force variations AF st ick = F m ax(i) — 
F m in(i) during a stick and AF s i ip = F ma x(i) — F min (i + 1) 
during a slip of a height H — 2.7D (126 g) of a miscellany of 
glass beads in the PMMA cylinder for three driving velocities: 
a) V = 150 nm.s -1 ; b) V = 750 nm.s -1 ; c) V = 1500 nm.s -1 . 



A. Stick-slip characterization 



For low velocities (V < 100 nm.s -1 ), a stick-slip mo- 
tion is obtained with the same properties as observed 
with monodisperse steel beads in a duralumin or brass 
cylinder. It is a regular stick-slip motion characterized by 
a maximum force before slipping that increases strongly 
when velocity is decreased, and a minimum force after 
slipping which does not depend on velocity (Fig. 122b ) . A 
coefficient K fJ, s (t s tick) can be extracted from F max with 
equation © for low velocities (V < 100 nm.s -1 ). Again, 
it shows a logarithmic aging of friction at the walls (Fig. 
122b). However, the behavior changes for higher velocities 
just before the transition to steady-sliding. The stick-slip 
motion is then less and less regular as the velocity is in- 
creased. Furthermore we observe structures in the force 
signal. Figure |2"31 evidences cycles of force increase dur- 
ing which the maximum and minimum forces increase for 
several stick-slip events, as the force amplitude AF s u c k 
during a stick is higher than the force amplitude AF s u p 
during a slip. Then, at the end of a cycle, there is a big 
slip event and a new cycle starts. This phenomenon is 
similar to the one observed by Albert et al. Q for the 
pushing of a stick in a granular material, and to the one 
observed for aluminium beads by Kolb et al. m the 
same display as ours but in 2D . It was however much 
more important in that last study in 2D: the force F max 



would then increase from 2A*" after a big slip event to 
20 N before the next big slip event. 

We analyse this phenomenon by studying the distribu- 
tion of force amplitudes defined as AF st i C k = F max {i) — 
F mm {i) for the i th stick event, and AF sUp = F max (i) - 
F m in{i + 1) for the following slip event (Fig. We ob- 
serve the evolution of these distributions with increasing 
velocity. AF st i C k distribution remains centered around a 
mean value; however, this distribution gets larger when 
the velocity is increased. For lower velocities (150 nm.s -1 
on Fig. |2"4"R). stick-slip is regular, and AF s i ip distribu- 
tion is regular around same mean value as AF st i C k (with 
nevertheless a larger distribution). When velocity is in- 
creased, the distribution is still predominant around same 
mean value as AF st i C k but many smaller values start to 
merge in the distribution (V = 750 nm.s -1 on Fig. 124b): 
these small values become as frequent as the larger ones 
for higher velocities (V = 1500 nm.s -1 on Fig. 124b). The 
distribution also enlarges to the higher values which cor- 
responds to occurrence of the big slip events. Note that in 
these experiments, the transition to steady-sliding regime 
occurs for velocity a velocity around V — 2000 nm.s -1 . 
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B. Analysis 

Several explanations can be proposed to account for 
this complex dynamics exhibiting a large distribution of 
energy release. In our point of view, the principal in- 
gredient for an explanation seems to be the progressive 
structuration of the packing between two big slip events. 
In the Janssen model picture, the progressive increase 
of F max in a cycle corresponds to progressive increase of 
Janssen's constant K. The K value increases from an ini- 
tial value Kq at the beginning of a cycle up to a maximum 
value K ma x at the end of the cycle, and is reinitialised 
by a big slip event to value Kq. Let us recall that K 
was found to increase with density in another report |15| . 
A consequence of this dependence was shown in section 
IIII A II Therefore, a way to obtain the phenomenology is 
to start from a loose packing at the beginning of a cycle; 
if each small slip event densities the packing (as can be 
expected from compaction under vibration), it leads to 
a small increase in K , so that the minimum force F m in 
at the end of the slip is higher than the force F m i n at 
the end of the preceding slip; the maximum force F max 
at the end of the stick is also higher than the force F max 
at the end of the preceding stick. Sometimes, big events 
occur; they probably correspond to convective motion at 
the walls, as observed by Kolb et al. in 2D ^lj, which 
loosens the packing, which then gets the initial density 
of the cycle. 

For the interpretation we give, the stick events last a 
constant time as the aging time of contacts at the walls 
is still controlled by the velocity. This explains why the 
distribution of AF st i C k is still centered on mean value 
contrary to AF s n p . 

But still a criterion for big slip events has yet to be 
found; it may be the level of stresses, or a dilatance effect 
at high density causing shear band and convection rolls. 
The regularity of stick-slip motion at low velocities can 
be understood if the criterion for big slip events is at an 
upper stress level. A slow driving velocity leads to long 
aging of contacts at the walls so that the maximum force 
Fmax is always greater than the threshold to obtain a big 
slip. All slip events are then big slip events in this case. 

This transition from a simple steady stick-slip to a 
complex dynamical regime with a large distribution of 
elastic energy release is an interesting issue and we plan 
to pursue further the investigation trying to clarify in 
particular the role of disorder and polydispersity in the 
phenomenology. 

V. CONCLUSION 

We presented an experimental report on the dynami- 
cal behavior of a granular column pushed vertically. We 
first investigated the case of a monodisperse assembly 



of steel beads and we observed at low driving velocities 
(1 < V < 100 fim/s) a steady behavior such that the 
pushing forces increase roughly logarithmically with ve- 
locity. On the other hand, at very low driving velocity 
(V < 1 iim/s), we evidenced a discontinuous and hys- 
teretic transition to a stick-slip regime characterized by a 
strong divergence of the maximal blockage force when the 
velocity goes to zero. All this phenomenology is strongly 
influenced by surrounding humidity: generally, higher 
humidity level increases strongly the resistance to push- 
ing. Finite elements numerical simulations were used to 
confront experimental results to a modelling of the gran- 
ular packing as an isotropic elastic medium. Then, we 
showed that a simple Janssen's model is a fair base for 
analysis as it provides the correct physical interpreta- 
tion for the pushing resistance. This Janssen's model 
was used to extract an important mechanical parameter 
combining the effects of stress redirection and wall fric- 
tion but there is an inherent difficulty to isolate clearly 
the various contributions either coming from bulk reor- 
ganization or from the surface friction properties. 

Using different column materials and measuring di- 
rectly the friction of a grain with the wall, we accumu- 
lated several evidences leading us to conclude that the 
force dependence with driving velocity and humidity is 
strongly related to the bead-wall friction properties: (i) 
in the steady limit grain/wall tribology measurements 
show a friction force increasing with humidity and ve- 
locity, (ii) in the stick-slip regime, the blockage enhance- 
ment can be related to humidity induced aging of the 
bead-wall friction, (iii) the hysteretic transition mecha- 
nism from stick-slip to steady-sliding is similar to the one 
observed in solid-on-solid experiments |27J, |28j . 

In spite of a dominant surface effect, we could also iden- 
tify contributions of bulk structurations. For example, 
we related the transitory part of the response to pushing 
to a dependance of the coefficient of redirection between 
horizontal and vertical stress with packing fraction. Also 
a clear difference of the mechanical parameters extracted 
from pushing and pulling experiments shows a contribu- 
tion from restructuration of about 20%. 

A second system, made of polydisperse assemblies of 
glass beads, was investigated. We emphasize the onset 
of a new complex dynamical behavior, i.e. the large dis- 
tribution of blockage forces evidenced in the stick-slip 
regime close to the transition. 
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